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SUMMARY 


The "buckling of a simply supported rectangular flat plate 
under combinations of shear and direct stress was investigated by 
means of an energy method. The critical combinations of stress for 
severe], length-width ratios were determined to an accuracy of about 
1 percent by the use of tenth-order determinants in conjunction with 
a modified matrix iteration method. Curves were drawn which can be 
used to obtain the critical stress combinations for the case of 
interaction of shear and longitudinal direct stress and for the case 
of interaction of shear and transverse direct stress. 


INTRODUCTION 


The problem of buckling of plates under the action of more than 
one stress has been given considerable attention. For plates 
subjected to two stresses the critical comb inations are usually given 
by means of interaction curves, that is, curves that can be used to 
obtain the value of one stresB required to produce buckling when a 
given value of another stress is also present. The stresses are 
ordinarily given nondimens ionally in terms of either stress ratios 
or stress coefficients. 

In reference 1 the Interaction curve for infinitely long flat 
plates under combined shear and longitudinal direct stress was 
shown to be very nearly a parabola as indicated in figure 1 of the 
present paper. In this figure R g is the ratio of shear stress 

present to the critical stress In pure shear and. R x is the ratio 

of the longitudinal direct stress present to the critical, stress 
in pure longitudinal compression • (All symbols are defined in 
appendix A.) This curve Ib shown in figure 1 of reference 1 to 
a T Ply to plates having edges either simply supported, .clamped, or 
elastically restrained against, rotation. Reference 2 demonstrates 
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■'hat- the parabola also applies ; to a high degree of accuracy, to 
nearly square plates with simply supported edges- On the other 
hand, for infinitely long plates loaded in shear and transverse 
direct stress, the interaction carve assumes a different form, 
which is shown for simply supported edges in figure 2 taken from 
figure 3, reference 3* The form of the interaction curve for shear 
and transverse direct stress appears, therefore, to change markedly 
as the length-width ratio of the plate increases from 1 to c©« 

The purpose of the present paper is to investigate this transition 
and, in addition, to determine whether any appreciable change In 
the form of the interaction curve for shear and longitudinal stress 
occurs as the length-width ratio of the plate increases from 1 
to oo. 


RESULTS AND DISCUSSICtT 


The results discussed herein are based on the theoretical 
solution presented in appendix B. The numerical results, computed 
by means of the matrix iteration method described in appendix C, 
are believed to be accurate to within about 1 percent. Interaction 
curves are presented for the buckling of simply supported plates 
having length-width ratios in the range of 1 to 4 and subjected to 
shear and longitudinal stress end to shear and transverse stress. 


Shear and Longitudinal Stress 

The interaction curves for shear and longitudinal direct stress 
are shown in figure 3 in terms of stress ratios for simply supported 
flat rectangular plates having length-width ratios of 1, 2, and 4. 

One of the two curves given for each length-width ratio represents 
a buckle pattern that is symmetric about the center point of the 
plate. The other curve represents a buckle pattern that is anti- 
symmetric about the center point of the plate. Those parts of the 
two curves that are governing parts (the parts that give lower 
values of one stress for a given value of the other stress) are 
drawn with solid lines . The purpose of presenting two curves 
instead of only the governing parts of either curve is to indicate 
where the cusps occur and to show the abrupt change from a symmetric 
to an antisymmetric buckle pattern that results in these cusps. 

For purposes of comparison, points from the parabola represented 
tho following equation for the interaction curve for a length- 
width ratio of co 

R 2 + R = 1 (seo reference l) 

S X 
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are included in each graph in figure 3 . The curves for the 
rectangular plates do not differ visibly from the parabola; this 
close agreement indicates that for shear and longitudinal direct 
stress of simply supported rectangular plates the parabolic 
interaction formula is substantially correct. 


Shear and Transverse Stress 

The interaction curves in stress-ratio form for shear and 
transverse direct stress for simply supported rectangular plates 
having length-vidlh ratios of 2, 3, and 4 are shown in figure 4. 

The curve for a length -width ratio of 2 is nearly a parabola in 
the compression range but deviates considerably from the parabola 
in the tension range. The curves for the higher length -width, 
ratios deviate considerably from a parabola in both compression and 
tension and incline toward the curve for infinitely long plates- 
In figure 5 the curves of figure 4, together with curves for plates 
of length-width ratios of 1 and co, are plotted on the same 
graph. The vertical straight-line part of the interaction curve 
for a length -width ratio of co corresponds to buckling of the 
plate as an Euler strip. (See reference 3*) 

In figare 6 the transition in the form of the interaction 
curve as the length-width ratio of the plate changes from 1 to co 
is shown in terms of buckling stress coefficients. Tho assumption 
that a long plate is infinitely long is seen to lead to a conservative 
estimate of the buckling stresses; if the length-width ratio is 
greater than 4, only a small error is involved in this assumption. 

The combinations of stress coefficients that result in buckling 
and the corresponding deflection coefficients are given in table 1. 


CONCLUSIONS 

The results for the critical combinations of direct stress and 
shear of simply supported rectangular flat plates, computed by an 
energy method summarized in tables and graphs, show that: 

1. For shear and longitudinal direct suress the interaction 
curve for all length-width ratios investigated is substantially a 
parabola for which the equation in terms of stress ratios is 

' 
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where B is the ratio of shear stress present to the critical stress 

9 

in pure shear and is the ratio of longitudinal direct stress 

present to the critical stress in pure longitudinal compression. 

2 . For shear and transverse direct stress, the shape of the 
interaction curve depends on the length -width ratio of the plate. 

For square plates the interaction curve is very nearly a parabola, 
the equation for which is given in the preceding paragraph. The 
interaction curve for plates having a length-width ratio of 2 is 
nearly parabolic in the compression range but is close to the hnovn 
interaction curve for infinitely long plates in the tension range . 

In the range of length -width ratio from 2 to 4 the parabolic equation 
does not hold even in compression; therefore the curves given in 
the present paper should be used. At length-width ratios greater 
than 4, the interaction curve approximates the interaction curve for 
a length-width ratio of oo . 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., November 8, 1946 
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.APPENDIX A 
SYMBOLS 


a 

b 

t 

w 

D 


length of plate 
width of plate; b S & 
thickness of plate 
deflection normal to plate 

flexural stiffness of plate 



x 

y 

B 

E 


longitudinal coordinate 
transverse coordinate 
symmetrical matrix 
Young’s modulus for material 


mn 


^mn 


mn 

Ci 


7 


T 

TL 


deflection-function coefficient 
element of matrix B 

diagonal terms in stability determinant 

mth unknown in set x-^x^, . . .2^ 
integers 

arbitrary coefficient 
longitudinal compressive stress 
transverse compressive stress 
shear stress 

longitudinal direct- stress ratio 
transverse direct- stress ratio 
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shear-stress ratio 

longitudinal compressive -stress coefficient 
transverse compressive-stress coefficient 
shear-stress coefficient 
length -width ratio (a/h) 

Poisson's ratio for material 

characteristic value of matrix B 

ith characteristic value of matrix B 

column matrix with elements x-^Xg,*-*^ 

modal column associated with ith characteristic value 

nth approximation to 2,- 

column matrix constructed so that each element is the mean of 
the corresponding elements of and ^ 


critical (used as subscript) 
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APPENDIX B 


THEORETICAL SOLUTION 
Method of Analysis 


The problem of the stability of finite rectangular flat plates 
under combined shear and direct stress 5s solved by use of the 
Rayleigh-RItz method. (See reference 4.) The deflection function 
Is expressed exactly by means of a two-dimensional Infinite 
Fourier series. The application of the Eayleigh-Ritz method results 
in an infinite set of homogeneous linear equations in the infinite 
nu-mber of unknown Fourier coefficients. Solutions of these equations 
that give Fourier coefficients not all equal to zero exist only for 
those combinations of shear and direct stress for which the buckled 
plate is in neutral equilibrium. 

The algebraic equations obtained herein are equivalent to 
the equations used in reference 2 , In order to obtain accurate 
results for relatively long plates, however, a method of solution 
of these equations different from the method of reference 2 was 
employed in the present paper. In reference 2, two sixth-order 
determinants (one determinant corresponding to a symmetrical and 
the other determinant to an anti symmetrical buckle pattern) were 
expanded and the resulting polynomials were solved for the critical 
stress combinations. Each solution involved the use of only 
six terms in the Fourier expansion of the deflection function. In 
the present paper 10 terms were used in the expansion of the deflection 
function. The corresponding set of 10 simultaneous equations was 
solved by the matrix interation method (references 4 and 5) modified 
in the manner described in appendix C. This method has the advantages 
that, whereas a very accurate solution requires a great deal of 
labor, a good approximate solution can be obtained without much 
effort, and the Fourier coefficients as well as the buckling loads 
can be evaluated. A knowledge of the Fourier coefficients for a 
given loading condition is used to determine whether the best 
choice of equations va3 made and also to assist in the choice of the 
10 most Important equations for a similar loading condition. 


Accuracy of Results 

Difficultly is usually encountered in evaluating the discrepancy 
between the buckling load determined by an exact, solution of the 10 
most important equations and the true value of the buckling load. 
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On the assumption that the use of the 10 most Important simultaneous 
equations of the infinite set to find the critical shear stress for 
a given compression stress vould in all cases result in an error of 
not more than a few tenths of 1 percent, the iteration (see appendix C) 
was carried out until the error in the solution of the equations used 
was estimated to "be less than 0.5 percent. The shear-stress 
coefficients presented in table 1 of the present paper therefore 
are believed to be in error by not more than 1 percent. 


Solution 

The critical stresses are determined on the basis of the principle 
that during buckling the elastic -strain energy stored in a structure 
is equal to the work done by the applied loads . For the case of a 
rectangular flat plate under loads applied in the plane of the plate, 
this equality becomes, when the coordinate system is that of figure 1 , 



(equivalent to equation 210, reference 6) 


Equation (Bl) can be rewritten in terms of nondimens icnal stress 
coefficients as follows: 
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where 


k 1§S 

x x « 2 d 


k .8 s|i 
• 7 7 « 2 b 


k = T 


felt 

n 2 D 


The procedure used in solving equation (332) involves substituting 
for w a function of x and y that satisfies the edge conditions 
and can he adjusted so as to closely approximate the buckling con - 
figuration. For any case in which the value of w 'is 0 at all the 
edges, the term with the coefficient -2(l - g) can be shown -to 
vanish (reference 7) . A series of terms with arbitrary coefficients 
is used to represent w, and the coefficients are determined by 
the Rayleigh-Bitz method. A general form for the deflection w is 


CO CO 

_ 1- a sin sin -~£ (B3) 

ro.ss.1 11=1 a b 


Equation (B2) is solved for the case of buckling under shear and 
transverse stress and for the case of buckling under shear and 
longitudinal stress, 

Shear an d trans v erse direct stress If the value of is 

set equal to 0 and the express ion for w Is substituted in' 
equation (B2), the following equation is obtained: 
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(E4) 

vhere m t p and n i q are odd numbers . 

The coefficient a must "be chosen to make the value of k 
mn 3 

a minimum. This procedure results in the set of homogeneous 

linear equations represented by 


ran 


TU 


2 p 2\ S 

n 13 i 


. 2J+ 

V p _ 




n2 


i v 

p=l q=l 


ronpq 


(m? - p2)(n' 


? . 


q 2 ) 

(B5) 


0 


where 

»-6 

m = 1, 2, 3, • • • 

n = 1, 2, 3, . . . 

and m t p and n ± q are odd numbers . 

This set of equations may be divided into two groups which are 
independent of each other, one group in which m i n is odd 
(antisymmetric buckling) , and one group in which m * n is even 
(symmetric buckling) . Ten equations in 10 unknowns were solved for 
each group by the iteration method explained in appendix C. 
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A representative determinant in terms of the coefficients for 
the group of equations in which m ± n is even is 




a ll 

a l3 

a 22 

1 — 1 
O'! 

a 

a 15 

a 24 

a 33 

a 42 

a 5l 

a 44 

m=l. 

n^l 

L 11 

0 

4 

0 

0 

8 

0 

8 

0 

16 




9 



45 


45 


225 

m=l, 

n=3 

0 

L U 

-4 

0 

0 

8 

0 

. 8 

0 

16 




5 



7 


25 


35 

m=2, 

n=2 

4 

4 

■ L 22 

4 

20 

0 

36 

0 

20 

0 



9 

'5 

5 

63 


25 


63 


m=3. 

n=l 

i 0 

• 0 

_4 

L 3i 

0 

8 

0 

8 

0 

16 





* R 
V 


25 


7 


35 

m-1, 

n=5 

0 

0 

20 

0 

L,- 

40 

0 

8 

0 

16 





63 


1? 

27 


"63 


"27 

m=2 } 

n«4 

8 

8 

0 

8 

40 

l 24 

_71 

0 

8 

0 




7 


~25 

"27 

35 


63 


di=3. 

n=3 

0 

0 

3.6 

0 

0 

.71 

L 33 

-72 

0 

144 





25 



35 

35 


49 

m=4, 

n=2 

8 

45 

8 

’25 

0 

8 

7 

8 

"63 

0 

_72 

35 

L 42 

4o 

'27 

0 

m=5. 

n=l 

0 

0 

20 

0 

0 

8 

0 

4o 

L 

0 





'63' 



'63 


'27 

51 


n=4 

16 

16 

0 

16 

16 

0 

144 

0 

0 

L 44 



225 

39 


35 

'27 


"49 




where 


inn 


32k s p3 


(m* + n% 2 ) 2 - yV 
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A typical determinant in terms of the coefficients for the 


group of equations 

in which m 

+ n is 

odd 

is 





a i2 

S 21 

a l4 

®23 

a 32 

e 4l 

a 25 

a 34 

a 43 

8 52 

m=l, n=2 

L 12 

4 

9 

0 

4 

5 

0 

8 

45 

20 

63 

0 

8 

25 

0 


ra=2, n=l 

4 

9 

L 21 

3 

45 

0 

4 

5 

0 

0 

8 

25 

0 

20 

63 


m-1, n=4 

0 

8 

L 14 

3 

0 

16 

4o 

0 

16 

0 




’45 

7 


"225 

21 


"35 



th-2, n=3 

4 

0 

8 

L 23 

.3.6 

0 

0 

72 

0 

4 



5 


" 7 

25 



35 


7 


m-3, n=2 

0 

4 

0 

.36 

l 3 2 

3 

4 

0 

22 

0 

(B7) 



5 


25 

’7 

"7 


35 



m=4, n=l 

3 

0 

l6 

0 

8 

L 4i 

0 

16 

0 

40 



"45 


“225 


"f 


'35 


27 


m=2, n=5 

20 

0 

4o 

0 

4 

0 

^5 

8 

0 

100 



63 


21 


7 


"3 


441 


m=3, n=4 

0 

3 

0 

72 

0 

16 

8 

L 34 

144 

0 




25 


35 


35 

" 3 

49 



m=4, n=3 

8 

25 

0 

16 

35 

0 

72 

35 

0 

0 

144 

49 

L 43 

8 

’3 


m=5, n=2 

! o 

! 

\ 

i 

20 

63 

0 

4 

’ 7 

0 

4o 

27 

100 

441 

0 

8 
" 3 

L 52 


In general, the method is 

to choose numerical 

values of 

p 

and k , 
y’ 

set 

each 

determinant 

equal to 0, 

and S' 

olve for the 

lowest 

value of 

V 

The 

lower 

of the two values 

of k s 

found 

from 

the 


two determinants establishes the critical shear stress for a panel 
with length-width ratio p under the transverse stress given by k y . 
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Shear an d longitudinal direct stre ss If k Is set equal 
to 0 in equation (B2) and the same procedure as^tho procedure 
for the combination of shear and transverse stress is carried out, 
the set of homogeneous linear equations given by the Rayleigh-P.itz 
method Is represented by 


mn 


| (rr? + n 2 p 2 ^ 


o O 

k m *fr 
x p 


+ 32kg£: 

«2 


\ 

z_.. 


N 2- 
Z_ a 


p=i q=i 


pq 


(m 1 


,2 - 


mn pq 

p 2 )(n 2 - q 2 ) 

(Bfl) 


vhere 
m = 1,2, 
n = 1, 3 > • * • 

and m + p and m + q are odd numbers . 

The determinants set up from these equations are the same as the 
determinants (b 6) and (B7) for symmetric and aixtisyimaotric buckling, 
respectively, except that the diagonal terms are 


■^mn 


32k s p3 


i ^m 2 + n 2 ^ 2 ^ 2 - k^m 2 



In general, the method is to choose numerical values of (3 
and k x , set each determinant equal to 0, and solve for the 

lowest value of k a . As in the determination of the shear and 
transverse direct stress, the lowest value of k Q establishes the 

critical shear stress for a panel with length- width ratio P under 
the longitudinal stress given by k^. 
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APPENDIX C 

MATRIX ITERATION METHOD 


If the matrix iteration method is used in attempting to solve 
the set of simultaneous linear algebraic equations associated with 
the buckling of flat plates where shear is present, the conventional 
iteration process as described in reference 4 does not converge . 

The reason for this nonconvergence, the modification in the iteration 
process used in the present paper to obtain convergence, and the 
method of choosing the best finite set of algebraic equations to use 
are described in the following paragraphs . 


Conventional Iteration Process 

The matrix iteration method is described in the present paper 
in terms of the following simplified symmetrical equation: 


CO 

Z ^rnn y n = ^m 

n=l 


(Cl) 


where 

m = 1 , 3 , ■ • ■ 

and 

b = b 

mn nm 

This set of equations is equivalent to the set of equations (B5) 
from which the critical shear stress is determined if 


m 


- 




+ , 1 2 3 2 ) 2 


- kyiV 


(C2) 


(where m is a different integer for each different combination 
of i and J ) , 


mn 



ijpq 

(i 2 - p 2 )(j 2 " q 2 ) | (i 2 + J 2 P 2 )‘ 


- v Sp f] 1/2 


(pq) 


1/2 

(C3) 
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(where n hears the same relation to p, q as m hears to 1, j), 
and 

x = -I- (C4) 

k s 

Equations (Cl) can similarly he shown to he a generalization of 
equations (b3) . 

The problem is to find the highest value of X which permits 
nonvanish ing solutions to he found for equation (Cl) . In order to 
avoid the use of too many suffixes, matrix notation is used, so that 
equation (Cl) becomes 

B£=X2, (C5) 

where B is a square matrix and 2, is a column vector. 

The matrix iteration method (reference f ) consists in taking 
some arbitrarily assumed values for the set of values for § 
and in calculating the left-hand side of equation (C5) in order to 
obtain an imnroved set of values for (Division by X is 

unnecessary because the solution of homogeneous equations can be 
determined only to within an arbitrary multiplicative constant.) 

These new values are reinserted in the left-hand side of equation (C5) 
to obtain further improvement in the values, and the process is 
continued until the ratios of the components of 2, are not 
appreciably changed by further Iterations. The value of X is then 
given as the ratio of the last value found for £ to the preceding 
value . 


The basis for this method can be seen from the following 
discussion. According 'to matrix theory, if tho matrix B is 
of Nth order, there are N values of X satisfying equation (C5) • 
Let these values he called X^_, X2 , . ••, Xjj with the order so 

arranged that | X-jJ > | Xg J > ••• > |Xjj| . For each value of X^ there 
is a corresponding solution for x, which may be called 2,^. Then 


= X& 


(C 6) 
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Ary initial assumed solution £^0) can "be expressed in terms 
of the N true solutions ?.jj as follows: 


C i S i (C7) 

i=l 

where the values of the coefficients are in general not known. 

The result of one iteration, obtained by using equations (C6) 
aM (C7), is seen to be 


JD 

S 



N 

= 

i«l 


N 


\ 

/ 

i=l 


c i x A 


The result of n iterations is 



H n. 

T. c i*A 

i»l 


The next step is to facftor out the nth power of X-^, that is, 
the largest. value of X: 


(n) 

£ 



ic^i 



(C8) 
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Equation (c8) shows that as n increases becomes more 

and more nearly a constant times the rate of convergence 

\ 

depending on the smallness of the ratio t— (i = 2, 3, • . • , N) . 

k l 

The largest value of X. can then be found as 

g(n) 

x i ■ r^u 

S 

Modification for Shear Buckling Problems 

In plate buckling problems in which shear i3 present, the 
critical shear stresses occur in pairs which are equal in magnitude 
but opposite in sign. For such a problem Xg = -Xq, and equation (C8) 
may be written 





(~D n c £ 2 + 




(C9) 


Equation (C9) implies that unless Cg is by chance equal to 0, 
no amount of iteration will result in convergence. 


A simple expedient can be used, however, to produce convergence. 
Since equation (C9) shows that as n increases g( n ' oscillates 
about the true solution (constant times £q), an improved approxi- 
mation can be obtained by constructing ^ such that each 
component is some kind of mean of the corresponding components of 

and Immediate convergence results from the use of the 

geometric mean (but not from the use of the arithmetic mean) after 
one iteration if only two equations of the set for a rectangular plate 
in shear are used 0 In the present paper, therefore, in winch 


10 equations were used, the geometric mean was employed except when 


Fourier coefficients of two successive iterations had opposite signs. 


in which case the arithmetic mean was used. 


The procedure adopted for obtaining the shear buckling stress 
was to use the matrix iteration method modified by taking the 
geometric mean after , every two iterations as described in the pre- 
ceding paragraph. In addition, by use of a method suggested in 
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reference 8, an Improved value for X]_ was obtained after the 
completion of n Iterations by means of the equation 


\ 


1 


F n) 5 (n ' X) 


Choice of Equations 

An exact solution of the critical shear stress for rectangular 
plates involves the use of an infinite set of equations in an 
infinite number of unknowns* Since attention must be confined to 
a finite number of equations - say, N - the ability to choose the 
beet N equations for the purpose is very desirable. 

A very useful (although not rigorously correct) guide to 
the best choice of the equations to be used may be obtained from a 
consideration of the accuracy of representation of the buckle 
deformation. The use of N equations in N unknowns implies that 
the deflection surface is being described in terms of N Fourier 
components, with the other components assumed equal to 0. 

The matrix iteration method yields the Fourier coefficients as 
well as the critical stress coefficient. The values found for these 
Fourier coefficients (where N was taken to be 10) were substituted 
in the following form: 


m + n even 


m + n odd 





& 17 

i 


a !8 




a 15 

e 26 



a l6 

a 27 



a 

13 

a 24 

a 

35 j 

1 ■ ■ _ 

a 

14 

a 

_!L 

a 

36 


*11 

.*22 

*33 

a 44 


a l2 

®23 

a 34 

a 45 



a 

31 

a 

42 

a 

53 


a 

21 

a 

32 

%3 

a 

54 




a 

51 

a 62 


a 4l _ 

a 52 

a 63 







~*6l_ 

V 




i 


*81 
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As a result of this substitution, values vere inserted in the 
10 squares corresponding to coefficients assumed not equal to 0, 
whereas no values vere substituted for the remaining squares. 

The absolute values of the coefficients vere observed to 
decrease in magnitude in a rather uniform manner as the distance 
from the largest coefficient increased. If a space in which no 
value was substituted occurs in one of the foregoing forms in a 
region where the neighboring computed coefficients are not small, 
appreciable error is usually incurred by the neglect of that term. 

In such cases, the buckling stress was recomputed with the coefficient 
in that space included and the smallest coefficient dropped from 
consideration. 
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TABLE 1 

9TRE3S COEFFICIENTS AND DEFLECT ION- FUNCTION COEFFICIENTS 





Shear and longitudinal 

direct stress 





s- 1 


Stress 




Deflectlon-funotlon coefficients a« i 

















k. 






1 ♦ j 

even 





Q 

i+J 

i.j 

a ll 

*13 

*22 

a 3i 

*15 

*24 

*33 

*42 

*51 

*35 


Even 

Odd 



L. 








4 

0 

7.05 

n 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3 

4.6a 

6.42 

B 

-.023 

.143 

-.017 

-.001 

.006 

.009 

.006 

-.001 

.001 

2 

6.62 

9.61 

B 

-.041 

.205 

-.035 

-.003 

.007 

.019 

.006 

-.002 

.003 

im 

6.11 


B 

-.057 

-253 

-.052 

-.004 

.006 

.028 

.005 

-.003 

.004 

1 ; 

9-35 

H.63 


-.070 

.293 

-.071 

-.005 

.005 

.036 

.003 

-.004 

.005 

-2 

11.56 


1 

-.091 

.360 

-.107 

-.006 

.002 

.057 

-.003 

-.006 

.007 

-4 

13.46 

14.76 

1 

-.107 

.**17 

-.145 

-.007* 

-.002 

.076 

-.011 

-.007 

,006 

* 




: 1 ♦ j 

odd . 1 




*12 

*21 

*14 

•23 

*32 

*41 

*25 

*34 

*43 

*5a 

4 

2. *9 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

3 

4.19 

3.29 

-.137 

1 

-.003 

-.016 

.149 

-.052 

-.001 

.005 

.007 

.014 


5.12 

4.66 

-.195 

1 

-.004 

-.035 

.214 

-.065 

-.002 

.006 

,015 

.014 

H 

5.93 

5.71 

-.239 

1 

-.004 

-.053 

.265 

-.109 

-.003 

.006 

.024 

.011 








1 ♦ J 

even 








*11 

*13 

*22 

*31 

V 

*24 

“33 

*42 

*51 

*35 

B 

6.59 

6.61 

1 

.032 

.343 

-.325 

-.002 

.005 

.046 

-.042 

-.010 

.004 

II 

7.69 


1 

.042 

.377 

-.266 

-.002 

. .005 

.056 

-.036 

-.014 

.006 

B 

9.04 

9.46 

1 

.052 

.407 

-.269 

-.003 

.003 

.069 

-.034 

-.016 

.006 









1 ♦ J 

odd 








a 12 

a 21 

*23 

*32 

*41 

*43 

*52 

•61 

*72 

*81 

4 

E9 

0 

0 

0 

0 

0 

X 

0 

0 

0 

0 

0 

3 

m 

2.9^ 

-.010 

-.224 

.004 

-.165 

1 

-.019 

.146 

-.066 

.025 

-.029 








i ♦ J 

even 








•n 

a 22 

‘31 

*33 

*42 

*51 

*53 

*62 

*71 

*91 

H 

4.10 

4.12 


.142 

1 

-.026 

-.2 52 

-.334 

.019 

.013 

-.026 

-.006 

B 


5.00 











0 

5.67 

5.77 


.'209 

1 

-.046 

-.324 

-.300 

.032 

,012 

-.032 

-.010 

-4 

7.98 

6.00 

-.315 

.316 

1 

-.090 

-.422 

-.292 

.057 

.025 

-.026 

-.010 


NATIONAL ADVISORY 
COMMITTEE FOR AERONAUTICS 














NACA TN No. 1223 


TABLE X 

STRESS COEFFICIENTS AND DEFLECTION-FUNCTION C0FFICIENT3 - Concluded 
3hear and transverse direct etreai 
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Fig. 3 
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Figure 3.- Interaction curves in terms of stress ratios R g and R^. 

for shear and longitudinal direct stress of simply supported 
rectangular flat plates having length -width ratios of 1, 2, and 4 and 
comparison with points from the curve representing the parabolic 


interaction formula 


R c 


R* - 1. 



; Rx - 




Fig. 4 


NAG A TN No. 1223 




A 

S 



Oy 


Antisymmetric buckling 
Symmetric buckling 


Ry 

NATIONAL ADVISORY x 
COMMITTEE FOR AERONAUTICS 


t 
l 
l 



i i ♦ t t t 


Figure 4.- Interaction curves in terms of stress ratios R s and R y 

for shear and transverse direct stress of simply supported 
rectangular flat plates having length -width ratios of 2, 3, and 4. 
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Figure 7.- Coordinate system for a rectangular flat plate. 
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